Theory of time-resolved optical spectroscopy on correlated electron systems 



OO 
O 

o 
< 



I 

C 

o 
o 



> 
in 
o 
o 

06 

o 

OO 

o 



X 



Martin Eckstein and Marcus Kollar 
Theoretical Physics III, Center for Electronic Correlations and Magnetism, 
Institute for Physics, University of Augsburg, 86135 Augsburg, Germany 
(Dated: August 7, 2008) 

The real-time dynamics of interacting electrons out of equilibrium contains detailed microscopic 
information about electronically correlated materials, which can be read out with time-resolved 
optical spectroscopy. The reflectivity that is typically measured in pump-probe experiments is 
related to the nonequilibrium optical conductivity. We show how to express this quantity in terms 
of real-time Green functions using dynamical mean-field theory. As an application we study the 
electrical response of the Falicov-Kimball model during the ultrafast buildup of the gapped phase 
at large interaction. 



I. INTRODUCTION 

Electronic correlations arc known to give rise to highly 
unusual phenomena such as heavy fermion behavior or 
the Mott metal-insulator transitions^ In recent years 
a new perspective for this field has been provided 
by various pump-probe spectroscopies, which can di- 
rectly track the time evolution of strongly interacting 
systems far from equilibrium. For example, the dy- 
namics of electrons in the vicinity of a Mott metal- 
insulator transition was investigated using time-resolved 
photoemission spectroscopy 2 - and time-resolved optical 
spectroscopy,Mi£i&iZ In these experiments, the sample is 
first excited by an intense laser pulse (pump); a sec- 
ond pulse (probe), which comes at a controlled time- 
delay, is then used to characterize the transient state by 
means of photoemission or optical spectroscopy. Pump- 
probe experiments with femtosecond time-resolution are 
now commonly used for the investigation of dynamics in 
molecules,— metals,— and semiconductors.— Recent devel- 
opment of shorter and shorter pulses has pushed the lim- 
iting time-resolution below 10 fs for optical frequencies^ 
and into the attosecond regime for pulses in the extreme 
ultraviolet^ 

For solids it is often a subtle task to distinguish the 
contribution of various degrees of freedom to a specific 
phenomenon. The Mott transition is induced by the 
Coulomb repulsion between electrons, but can occur si- 
multaneously with a change of the lattice structure, ob- 
scuring the primary origin of the phase transition. In 
time-resolved experiments, however, different degrees of 
freedom can be identified if they evolve on different time 
scales^ In particular, the lattice usually reacts much 
slower than the electronic system. Many phenomena that 
are already visible at low time resolution can be explained 
in terms of a two-temperature model^ 2 ^ which assumes 
that the electronic system is in thermal equilibrium at 
any given time, but may have a different temperature 
than the lattice. 

On the other hand, pump-probe experiments allow for 
an investigation of the electronic real-time dynamics. For 
example, two-photon photoemission spectroscopy can 
monitor the ultrafast thermalization of a pumped elec- 



tron gas in metals within several 100 fs^^ In semicon- 
ducting GaAs, the buildup of a screened Coulomb inter- 
action in the electron-hole plasma created by the pho- 
toexcitation of electrons into the conduction band has 
been tracked using time-domain THz spectroscopy^ In 
particular the latter experiment probes the true quantum 
dynamics of the state, which can no longer be described 
by a simple rate equation but requires the full many- 
particle Hamiltoniani 16 i 17 It would be very interesting to 
measure the electronic dynamics in strongly interacting 
systems, which may dominate, e.g., the ultrafast buildup 
of intermediate metallic states across insulator-to-metal 
transitions^^ 7 - or the melting of correlation-induced 
long-range order after an external perturbation.— The 
goal of this paper is to set up the framework for a micro- 
scopic description of time-resolved optical measurements 
in such strongly correlated electron systems. For time- 
resolved photoemission spectroscopy, the microscopic de- 
scription was recently derived in Ref . Hsl 

The microscopic formalism of isolated quantum many- 
body systems out of equilibrium was given independently 
by Baym and Kadanoffr 9 - and Keldysh3£ in terms of real- 
time Green functions. It provides the starting point for 
a nonequilibrium perturbation theory ; 21 ' 22 which is how- 
ever bound to fail for strong interactions. Dynamical 
mean- field theory (DMFT) ^ which becomes exact in the 
limit of infinite spatial dimension^ also applies to the 
non-perturbativc regime. DMFT self-consistcntly maps 
a lattice model onto an auxiliary single-site problem. The 
equilibrium theory has been instrumental in understand- 
ing many correlation-induced phenomena, such as the 
Mott transition, both for simple model systems^ and 
for real materials! 25 ' 26 Recently, DMFT for nonequilib- 
rium has been formulated in the framework of Keldysh 
theory^ It has been used to investigate the Falicov- 
Kimball mode l 28,29 under the influence of strong elec- 
trical fields ; 27 ' 30 as well as its relaxation over the metal- 
insulator transition after a sudden change of the interac- 
tion parameter i^i Similar investigations for the Hubbard 
model still require new techniques for the solution of the 
effective single-site problem. However, promising candi- 
dates for this task have been developed during the last 



years 
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The main purpose of this paper is to discuss the probe 
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process in optical spectroscopy in terms of linear re- 
sponse of a nonequilibrium state to an electromagnetic 
field. For this state, which might originate from the 
application of a pump pulse, the time evolution is as- 
sumed to be known from DMFT. The response is given 
by the two-time optical conductivity a(t, t'), that relates 
the current at time t to electrical fields in the sample 
at earlier times t' ^ For systems in equilibrium DMFT 
has already been successfully used to understand op- 
tical spectroscopy in correlated materials^ The stan- 
dard expression for the frequency-dependent conductiv- 
ity ct(cj) in DMFT— is quite simple and contains only 
single-particle Green functions, because vertex correc- 
tions to the current-current correlation function vanish 
for isotropic systems i 36 ^ 37 In this paper we derive an ex- 
pression for the two-time conductivity a(t, t') from non- 
equilibrium DMFT, which turns out to be a direct gener- 
alization of the equilibrium expression^ to Keldysh lan- 
guage. In particular, our derivation shows when the in- 
clusion of vertex corrections becomes mandatory in non- 
equilibrium situations, and under which conditions simi- 
lar simplification occur for a(t,t') as for ct(w). 

We then apply the theory to a simple lattice model for 
interacting electrons in a single band, 



H 



ija 



ij C ia C 3<? 



(1) 



Here cfj are annihilation (creation) operators for two 
species of fermions (cr = J,, "J") on lattice site i, which in- 
teract via a local Coulomb repulsion U. The first term in 
(fl]) is a tight-binding description of the electronic band. 
Eq. (fl]) is the Hamiltonian of the defines the Hubbard 
model if the hopping V? does not depend on the flavor 
cr, or the Falicov-Kimball model^ if one particle species 
is immobile (Vy = 0). Both models have a rich phase di- 
agram as a function of interaction and filling, including 
metallic, insulation and ordered phases. In the presence 
of electromagnetic fields [with scalar and vector poten- 
tial $>(r,t) and A(r,t)], the hopping amplitudes acquire 
Pcicrls phase factor o 38 i 39 



Vij = V,j exp 



— JdrA{r,t) 



(2) 



and a potential term — $(l?j,i)cj cr Cj (T is added to 
the Hamiltonian, where — e is the charge of an electron. 
Here and throughout a tilde indicates that the quantity 
is taken in zero external field. 

Nonequilibrium DMFT can potentially model the full 
pump-probe process by including the pump field explic- 
itly in Eq. @ . In the application of the general result to 
the Falicov-Kimball model we use an idealized nonequi- 
librium situation instead, where the "pumping" is an in- 
stantaneous event; we only have to know the excited state 
after the pumping, which is taken as initial state for the 



subsequent time evolution. This permits an investiga- 
tion of the relaxation between the various phases. For 
instance, we can start from a metallic state and follow 
the relaxation in the insulating parameter regime of the 
Hamiltonian. Below we model this situation by a sudden 
increase of the interaction parameter U. We therefore 
allow for arbitrary time dependence of all parameters U, 
\i and V? in the Hamiltonian {1}. 

This paper is outlined as follows. In Section [Til wc 
define the optical conductivity in nonequilibrium experi- 
ments, and discuss its relation to the reflectivity in time- 
resolved measurements. In Section [IIII we shortly review 
DMFT for nonequilibrium. We then derive the nonequi- 
librium optical conductivity in DMFT in Section HVl Fi- 
nally, in Section [V] we apply the theory to the Falicov- 
Kimball model and investigate the response of the system 
during the ultrafast buildup of the gapped phase at large 
interaction. 



II. TIME-RESOLVED OPTICAL 
SPECTROSCOPY 

To understand the results of time-resolved optical spec- 
troscopy it is necessary to know how weak electromag- 
netic pulses of finite length propagate through the sam- 
ple, which is not in equilibrium due to the applied pump 
pulse i 34 ' 40 ! 41 The current Sj is the linear response in- 
duced by the probe field, 



Sja(r,t)= / dt , a al3 (t,t')SE l3 (r,t'), 



(3) 



which defines the optical conductivity <7 a p(t,t') for sam- 
ples that are not in equilibrium. (Here and throughout 
a and /3 arc cartesian components of the vectors, and 
repeated indices are summed over.) Note that only the 
response ((3|) is linear in the probe field 5Ep(r,t'), while 
arbitrarily strong electric pump fields might be acting 
on the sample. The wavelength in optical spectroscopy 
is typically much larger than the lattice spacing of the 
sample, so that the linear response relation (J3]) is essen- 
tially local in space. On the other hand it is not local 
in time, and unless there is a clear separation between 
the time scales that govern the electromagnetic response 
and the relaxation of the nonequilibrium state, a(t,t') 
depends not only on the difference of its time arguments 
but on both t and t' separately. Of course a(t,t') is al- 
ways causal, i.e., it vanishes for t < t' . 

Knowledge of a(t, t') is sufficient to calculate the re- 
flected and transmitted pulses from Maxwell's equations, 
assuming that the induced current inside the sample is 
given by Eq. Q , 34,40,41 jj oweverj the relation to mea- 
surable quantities is more complicated than for samples 
that are in equilibrium. To illustrate we this consider 
a typical time-resolved reflection experiment, performed 
at normal incidence, on a sample that is infinite in the 
y-z plane (cf. Fig. [JJ). Outside the sample light propa- 
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FIG. 1: Time-resolved reflection experiment. For f — > — oo a 
probe pulse Eo(t,x) = zEo(t — td — x/c) propagates in +x 
direction without dispersion (upper panel). The sample is hit 
at times near t — td, and emits the reflected pulse E Te a(t, x) 
— zE Te &(t — td + x/c), which propagates in — x direction after 
leaving the sample (lower panel). 



gates without dispersion, so that we may write Eo(t,x) 
= zE (t-t d - x/c) and E rcfi (t,x) = zE Icfl (t-t d + x/c) 
for incident and reflected pulses at x — ► — oo, respectively. 
The functions Eq(t) and i? re fl(r) are centered around r = 
0, and t d is the probe delay. For simplicity we assumed 
cubic symmetry, such that the polarization direction z 
for both pulses is the same. We then define a generalized 
reflection coefficient r(t,t'),— 



E ic r{t) = J ds r(t d + T,t d + r - s)E {t 



s), (4) 



providing a linear relation between the two pulses. The 
full two-dimensional time-dependence of r(t,t') can be 
deduced from experiment by suitably choosing the pulse, 
and measuring at all possible pump-probe delays t d . 
However, if the optical conductivity a(t, t') depends on t 
and t' separately, then there is no simple relation to the 
reflection coefficient r(t,t')A^ This is evident from the 
definition Q , which shows that a sample which is not in 
equilibrium can modulate the pulse frequency. 

From now on we use an approximate form for r(t,t'), 
which is valid for reflection from a very thin slab (with 
thickness L — ► 0), such that the phase lag between the 
borders is negligible. In this case Maxwell's equations 
are easily solved, yielding^ 



r(t,t') = -a(t,t'). 
c 



(5) 



A more realistic description, which takes the finite thick- 
ness of the sample and its inhomogeneous excited state 
into account, requires the numerical simulation of the 
pulse propagation^ and of the inverse problem^ of ob- 
taining a(t,t') from r(t,t'). However, the treatment of 
such effects is beyond the scope of this paper, the goal of 
which is to calculate the optical conductivity a(t,t') mi- 
croscopically for an interacting many-body system that 
is not in equilibrium. 



III. DMFT FOR NONEQUILIBRIUM 

DMFT for noncquilibrium usually starts from ther- 
mal equilibrium at some early time t = t m i n i 27 i 31 For 
t > tmin the system evolves according to the Hamilto- 
nian ([1]), driven out of equilibrium if the Hamiltonian 
changes with time. Thermodynamic variables and op- 
tical response functions are obtained from the retarded, 
advanced and lesser real-time Green functions, 



^{t~t'){{c la {t),c) a {t')}) (6a) 
Gf ja (t,t>) =ie(t'-t){{c ia (t),cl(t')}) (6b) 
Gf ja (t,t')=i(cl(t')c i(T (t)). (6c) 

(Although retarded and advanced Green functions are 
in fact related by symmetry, both are given here for 
later reference.) The average (•) = Tr[p -] in Eq. © 
is over initial states at t — t m ; n , distributed ac- 
cording to the grand-canonical density matrix p oc 
exp[— /3H(t m i n )] at inverse temperature /?. The op- 
erators Ci a (t) = U(t, t m in)ci (T J7(tmin, t) are in Heisen- 
berg representation with respect to the full time evolu- 
tion U(t,t') = T f exp[-i It'diH(t)}. Using the Keldysh 
formalis m 20 ' 21 i 22 the Green functions ([6]) are then calcu- 
lated in terms of a more general contour-ordered Green 
function Gij, a (t,t') = — i(Tcc jer (t)ct.(t')} with time ar- 
guments on the contour C that runs from i m ; n to some 
larger time t m ax on the real axis, then from t max to £ m in, 
and finally to t m ; n — i/3 on the imaginary time axis. For 
the retarded, advanced and lesser components one has^i 



Gf j<7 (t,t')=Q(t-t')[GT j t(t,t')-G+-(t,t')] (7a) 
= G±t(t,t')-G+-(t,t') (7b) 
G^(M') = 6(t' -t)[G±-{t,t') GT+{t,t')] (7c) 

(7d) 
(7e) 



Gt-(t,t')-G--(t,t') 



Gfjait, t') — Gfj a (t, t'), 



where the superscripts ± indicate whether the first and 
second time arguments are on the upper or lower real- 
time branch of the contour, respectively. 

From now on we only consider translationally invari- 
ant nonequilibrium states, i.e., we assume that the Green 
function Gij ia (t,t') depends only on the difference Ri — 
Rj, with diagonal Fourier transform Gkcr(t,t'). This as- 
sumes that the electromagnetic fields do not depend on 
explicitly on position either, which is justified for experi- 
ments at optical frequencies, as discussed above in Sec.lHl 
We use a gauge with zero scalar potential $, for which 
electrical field is given by E(t) = —dtA(t)/c. The hop- 
ping amplitude Vij [Eq. ^] then also depends only on 
the distance Ri — Rj] its Fourier transform ek&{t) is given 
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£ka(t) = ^2 V i0 ex V>[ ik ( R j - #t)] = ^k+^A(t),a , (8a) 
3 

Ika = % eMiHRj - Ri)} , (8b) 

3 

i.e., e kcr (t) is obtained from the zero-field dispersion e keT 
by a time-dependent shift in momentum. 

The interacting contour Green function satisfies the 
Dyson cquatio n 21 i 22 

[{Q k l - E feff ) * G kcr ](t,t') = S c (t,t% (9) 

where T, kcr (t,t') is the contour self-energy and Gka(t,t') 
is the noninteracting Green function, whose inverse 

Q^(t,t') = 5 c (t,t)[i% + (/v - e k „(t))/h] (10) 

can be written as a differential operator on the contour. 
Here (/ * g)(t,t') = J c dtf(t,t)g(t,t') is the convolution 
of two functions along the contour, 8c{t,t') is the con- 
tour delta function [defined by J c di f(t)5c(i,t) = f(t)}, 
and <9f denotes the contour derivative.— The unique so- 
lution of Eq. @ is determined by antiperiodic boundary 
conditions for the contour Green functions in both time 
arguments . 21 ' 22 

The DMFT self-energy is local in space, i.e., H kcr is 
independent of k for a translationally invariant system. 
This approximation becomes exact in the limit of infinite 
spatial dimensions, 24 both for the equilibrium self-energy 
and the Keldysh self-energy— In DMFT the local self- 
energy E CT (i, t') and the local Green function G CT (i,i'), 

G c {t,t') = G Ua (t,t') = l]TG felT (M'), (11) 

k 

(N is the number of lattice sites in the sample) are deter- 
mined from an auxiliary problem in which the degrees of 
freedom at a single lattice site i are coupled to some un- 
known environment. The latter must be determined self- 
consistcntly, by solving the auxiliary problem together 
with the Dyson equation © . As the precise form of the 
local problem in terms of its many-body action does not 
enter into the derivation of the electromagnetic response 
below, we refer to previous work for further detailsi 27 ' 31 

IV. OPTICAL CONDUCTIVITY IN DMFT 

The current operator for the Hamiltonian (TTJ) is 
definedi^iii by the relation j(r) = -c5H/6A(r). Us- 
ing Eq. © , we obtain the current in the long wave-length 
limit as 

m) = (vj ddr ^ r y qr ) o > ( i2a ) 

= ^E^w G ^)' ( 12b ) 

ka 



the current vertex is given by 

Vka(t) = ft~ 1 9fce fcCT (t) = fir 1 <9fce fc+ _ fej4(t)j(T . (12c) 

and V is the volume of the sample. Although the 
response to arbitrarily strong fields is described by 
DMFT—, here we are interested in the linear current re- 
sponse to a weak probe field. We define the susceptibility 

x a pfrt') = s(j<*(t))/M t ')- (13) 

In the chosen gauge with E(t) = —d t A(t)/c, the sus- 
ceptibility Xafl(t, t') is related to the optical conductivity 
a aP {t,t') [Eq. ©] by 

OO 

° a p{t,t') = -cJ dtx^t). (14) 

v 

The susceptibility (fl~3|) is related to the current-current 
correlation function, which can be evaluated in analogy 
to the equilibrium case^ Here we prefer to take the 
derivative of (|12b|) directly, where the vector potential en- 
ters both in the vertex v k(T (t) and in the Green function 
G kr7 (t, t). This yields the diamagnctic and paramagnetic 
contributions to the susceptibility, 

Xap(t, = x^ a (i, + Xl™(t, , (15a) 

fcer ^ > 

*5UO = ^£^0^|M. (150 

ka ^ 

The paramagnetic contribution can be found from a vari- 
ation of the lattice Dyson equation ([5]), 

SG ha = -G ka * [6g£ - * G ka . (16) 

Some simplifications occur in the absence of 
anisotropies. We note that the second term in (|16|). 
containing the fe-independent self-energy, does not con- 
tribute to the fc-sum in Eq. (|15cp if, under inversion of k, 
(i) G kl j is symmetric and (ii) the vertex v ka is antisym- 
metric. These conditions are met by an isotropic system 
without external fields, and are therefore generally valid 
for systems with inversion symmetry in equilibrium. 37 
However, the isotropy may be lost when an initially 
isotropic system is driven out of equilibrium, e.g., when a 
current is induced by the electrical pump field. Further- 
more, the vertex (|12c[) is no longer antisymmetric when 
an electrical field is present in addition to the probe field, 
i.e., when the paramagnetic susceptibility (|15cj) is evalu- 
ated at A 7^ 0. Experimentally these anisotropic effects 
in otherwise isotropic systems show up as a dependence of 
the signal on the relative polarization of pump and probe 
pulses. However, when the anisotropy is caused entirely 
by the pump pulse, the inversion symmetry of G ka * 8H a * 
G k<J can be restored by averaging over the pump pulse 
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polarization. Then this term again drops out in (|15c|) . 
provided that Vk a {t) is antisymmetric (i.e., A(t) = 0). 
In order to study such anisotropic effects, vertex correc- 
tions contained in 5T, a must be taken into account (even 
for cubic lattices), by solving a Bethe-Salpeter equation 
on the Keldysh contour, with the irreducible vertex func- 
tion 5Ti a /5G a from the auxiliary single-site problem as 
input. 

In the following we only consider the completely 
isotropic relaxation between homogeneous phases, such 
that the vertex corrections (5E CT can be disregarded. 
Eq. (|15cp is evaluated at zero field, so that only the first 
term 5F ka (t 1 ,t 2 ) = -[G k<J * 6{Q^) * G ka ]{t 1 ,t 2 ) con- 
tributes to 5Gka in Eq. (fl6|) . This corresponds to keep- 
ing only the elementary bubble diagram for the current- 
current correlation function^ The two convolutions in 
SFka(ti,t2) collapse to a single one because [<5£7fcJ](£, £') 
(x Sc(t,t'). In order to obtain SGi (t, £) we take t\ = £ 
and £2 = t on the upper and lower branch of the contour, 
respectively [cf. Eq. (J7J)]. The contour integral is then 
transformed into an integral along the real axis, 



6F k a(t+,t-) = ^ J div ka {t)SA(t) 

— OO 

x [G + k t (*. mt (*> t) G+- (£, t)G~ (£, £)] (17) 

from which the optical conductivity a a p(t, £') can be read 
off. From Eq. (JT]). together with the relations G ka {t,t') 
= -G^(£',£)* and GL(£,£') = G^(f,t)*, we finally 
obtain the paramagnetic susceptibility 



X %%(t, £') = -2 X oE v^vllmlGUt, t')G< c (t', £)], 

(18a) 



ko- 



wherc xo = e 2 /(Vhc) and v ka = dktk.a/fo- The diamag- 
netic contribution follows directly from (|12cjl : 



Xo 



^ a (£,£') = l -fd(t - 0£(9kA^*«r)G&(*.t) • (18b) 



Eqs. (fTi)) and (JTHJ) constitute our final DMFT expressions 
for the optical conductivity (provided that anisotropic 
effects are disregarded, as discussed above). 

The optical conductivity (fH|) can be written as 

a a p{t,t') = K°|(£,£') + D aP (t)]e(t - £') , (19) 

i.e., it splits into its regular part 



(20) 



contribution 

D a p{t)= Mm a af3 (£,£') 

t' — > — OO 



= a^(t)-c / dt X ^(t,t) 



(21a) 



(21b) 



which does not depend on the time difference at all. In 
the latter expression, erj|(£) = — dt' Xa0 (*> *') is 
the weight of the delta function in Eq. ()18bj) . A finite 
Drudc contribution D a /}(t) 7^ indicates perfect metal- 
lic behavior, because it gives rise to a delta function at 
zero frequency in the partially Fourier-transformed opti- 
cal conductivity 



a a p(t, u) = Jds e^ +40) V Q/3 (£, t - s) (22a) 




iD a p{t) 
ui + iO 



(22b) 



Note that Eqs. (fT5| and (|2"0"|) can be checked by inserting 
equilibrium Green functions 



G£ CT (£, £') = -iQ(t - £') / duj A^y^'-V , (23a) 
GtrM = * I dujAkA^)f^)e lul{t '- t) , (23b) 



with the spectral function Ak a {ui) = — I m [G^ CT (w + 
«0)]/-7r and the Fermi function f(uj) = l/(l + e' 9w ), which 
depend only on time differences, into Eq. (|22|). Then the 
well-known expression for the regular part of the optical 
conductivity in equilibrium^ 



Re (uj) = ttcxo £ v^v^ 



which vanishes in the limit £' — > —00, and the Drude 



is recovered. 



V. PUMP-PROBE SPECTROSCOPY ON THE 
FALICOV-KIMBALL MODEL 

A. The Falicov-Kimball model in nonequilibrium 

In the remaining part of this paper we focus on a spe- 
cific electronic model, the Falicov-Kimball model. This 
lattice model describes itinerant (J,) electrons and im- 
mobile (|) electrons that interact via a repulsive local 
interaction The Hamiltonian is given by Eq. (p} 

with V}, = 0. The Falicov-Kimball model has been 
an important benchmark for the development of DMFT 
in equilibrium, because the effective single-site problem 
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for the mobile particles is quadratic and can be solved 
exactly^ It currently plays a similar role for nonequi- 
librium DMFT . 27 ' 30 ' 31 in particular since no appropriate 
real-time impurity solver is yet available for the Hub- 
bard model. In spite of its apparent simplicity the 
Falicov-Kimball model has a rich phase diagram contain- 
ing metallic, insulating, and charge-ordered phases^ In 
the following we fix the filling of both particle species 
(ni = n-\ = 1/2), and consider only the homogeneous 
phase without symmetry breaking. This phase under- 
goes a metal-insulator transition at a critical interaction 
U = [7 e ) 45 i 46 i 47 from the gapless phase at U < U c to the 
gapped phase at U > U c . 

Below we assume that the system is prepared in ther- 
mal equilibrium for times t < 0. Then the interaction 
parameter U is changed abruptly at t = 0. In this way 
we study the relaxation of the system in the insulating 
parameter regime, starting from a weakly correlated state 
(U < U c ). This mimics an experiment similar to the one 
described in Ref. [H, where the buildup of a weakly cor- 
related state is studied with time-resolved spectroscopy, 
starting from an uncorrelated state of electrons just after 
their excitation into an empty conduction band. Note 
that in this interpretation the state of the conduction 
band immediately after the pump pulse is the initial state 
for the relaxation process. 

The relaxation dynamics after such an interaction 
quench was recently investigated with DMFT using the 
exact Green functions Gki(t, t') of the mobile particles^ 
However, only thermodynamic observables were dis- 
cussed in Ref. [U with a special focus on their steady 
state value in the long-time limit. Here we consider 
instead hypothetical time-resolved experiments that are 
performed on the system during relaxation, i.e., we use 
the Green functions from Ref. [3l| to evaluate the opti- 
cal conductivity from Eq. (|18[) . Momentum summations 
in (|18[) are performed for a hypercubic lattice, taking 
the dispersion £k to be that of a semielliptic density of 
states^ p(e) = {2/nW 2 )VW 2 - e 2 (cf. App. EJ). The 
half-bandwidth W = 2 sets the energy scale, such that 
the critical interaction is U r = W = 2. 



B. Optical conductivity and reflected electrical 
field 



We study relaxation far in the insulating regime (U = 
6), starting from an initial metallic state (U = 1). The 
optical conductivity a(t, t — s) for this case is shown in 
Fig. [2k as a function of t and s. There are five regions 
[(i) to (v)] in this plot that we want to discuss in detail. 

In regions (i) [t < 0] and (ii) [t — > oo], a(t,t — s) de- 
pends only on the time-difference s, indicating that the 
system is in a stationary state. For (i) this is the initial 
equilibrium state, and for (ii) it corresponds to the final 
steady stated The Fourier transformation (|22a[) of the 
conductivity exhibits a broad peak at u> = 0, both for the 
initial state \&{t = 0,u>)] and the final state [a(t = oo,u>)) 
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FIG. 2: (a) Optical conductivity a(t,t — s) for the quench 
from the ground state at U — 1 (initial temperature T = 0) 
to U = 6 (n f =n c = 1/2, half-bandwidth W = 2). The unit 
of the conductivity is go = Na 2 e 2 W/(2h 2 V), where a is the 
lattice constant. In the region above the upper dashed line, 
t — s < 0. Below the lower dashed line the relaxation is es- 
sentially complete, (b) Fourier transform p2a[) of the optical 
conductivity in the initial and final stationary state, and for 
an equilibrium state at U — 6, with the same excitation en- 
ergy relative as the final state (T = 2.070). (c) Diamagnetic 
contribution (|18bl) to the susceptibility. 



(cf . Fig. (2b) • This clear indication of metallic behavior of 
the final state may seem surprising, since the interaction 
is far above the critical interaction U c . However, a finite 
DC conductivity should be expected because the final 
state is highly excited with respect to the ground state 
at U = 6. In fact, the excitation energy corresponds to 
an effective temperature T = 2.070, for which the equilib- 
rium DC conductivity er(0) is already quite sizable even 
at U — 6 (dotted curve in Fig. [^b). However, er(0) is 
still considerably lower compared to d{t = oo, 0). This is 
a signature of the incomplete relaxation in the Falicov- 
Kimball model: The system does not relax to thermal 
equilibrium, but reaches a non-thermal steady state, as 
shown in Ref. [3l| for thermodynamic quantities. In the 
present context we find that the electromagnetic response 
of the non-thermal final state combines some features of 
the insulating state (a peak around to = 6 due to excita- 
tions across the gap) with a sizable DC conductivity. Full 



7 



thcrmalization is expected only due to coupling to fur- 
ther degrees of freedom or further hopping or interaction 
terms that are not contained in ([1]). 

For t—s < and t > [region (hi) in Fig. [2^], a(t,t—s) 
determines the current after the pumping at t = caused 
by an electrical held applied to the sample before the 
pumping. It thus measures a combination of the electro- 
magnetic response of the initial state and the subsequent 
decay of the induced current for t > 0. By contrast, 
in region (iv) in Fig. [2K it describes the response of the 
noncquilibrium state alone, and hence gives direct insight 
into various relaxation processes. True noncquilibrium 
dynamics can be observed only when both t — s and t are 
smaller than some relaxation time r stat , after which the 
response is stationary, i.e., when a(t, t—s) depends on s 
only. In the present case the relaxation is virtually com- 
plete after only a few times of the inverse half-bandwidth 
( Tstat ~ 8/W = 4, below the lower dotted line in Fig.[2K)- 
Therefore the relaxation time and the time scales of the 
electromagnetic response, which is set by the decline of 
a(t, t — s) for s — * oo, apparently have the same order of 
magnitude. 

In spite of this very fast relaxation nontrivial tran- 
sient behavior can be observed before the stationary state 
is reached. Consider a(t, t — s) at s = 0, which tra- 
verses almost two damped oscillation cycles with an ap- 
proximate period 2'xh/U before reaching its final value 
(Fig. Recall that cr(t,t) is given by the delta func- 
tion weight <r dla (t) of the diamagnetic susceptibility (|18bp 
[cf. Eqs. ([20]) and ([21])] . These oscillations arc the hall- 
mark of dynamics that are dominated by a Hubbard-type 
density interaction such as [/^^ n^n^. In fact, when 
the Hamiltonian is given only by this interaction term, 
the time evolution-operator exp[itU '^2 i n^jj itself is 
time-periodic^ and oscillations should therefore be vis- 
ible in all nonlocal quantities. These so-called collapsc- 
and-revival oscillations were first observed and described 
in experiments with ultra-cold atomic gases,— where the 
Hamiltonian of the system can be designed in a controlled 
way. 

Finally we note that the conductivity a(t,t — s) van- 
ishes in the limit s — » oo , i.e., the Drude weight (|2~lj) 
vanishes for all times [region (v) in Fig. ■ This is well 
known for the Falicov-Kimball model in equilibrium: 29 
unlike in the Hubbard model^ the mobile particles do 
not form a perfect metal even at T = because of the 
disordered background of immobile particles. Mathemat- 
ically the vanishing of D a p(t) is due to the cancellation 
of the two terms in (|2"T|) . Since each of them has a non- 
trivial time dependence (cf. Fig. this cancellation 
represents a strong check for our numerical evaluation of 
the conductivity. 

To illustrate the relation of the optical conductivity 
to time-resolved THz experiments, we use the simple 
expression (|5|) for the reflection coefficient, and calcu- 
late the reflected field -Ercfl (t; id) according to the def- 
inition @, using a single cycle incident pulse Eq(t) — 
sin(r) exp(— 2r 2 ). The result is shown in Fig. [3] For short 




FIG. 3: Result of an idealized spectroscopy experiment 
(cf. Fig.[T]): (a) Incident pulse, (b) Reflected pulse E le a(r;td) 
[from Eqn. Q and (J5J] , for a delay td of the incident pulse 
with respect to the start of the relaxation at t = 0. The re- 
gion below the diagonal dotted line is not influenced by the 
quench at all. Above the horizontal dotted line (td > T s t a t = 
4) the reflected signal is converged. For td < T sta t, at least 
one revival peak at td — 2nh/U is clearly visible (crosses). 



delay times td between the incident pulse and the pump- 
event at t = 0, the profile of the reflected field depend 
strongly on t^. On the other hand, for times td > T stat , 
the relaxation is essentially complete, and -E re fl(T) has 
developed a longer oscillating tail. This general behavior 
is also seen in the experiment of Ref. [HI. In Fig. [3] the os- 
cillations in i? re fl(r) as a function of r are characteristic 
of the gap in the final state. Furthermore, the above- 
mentioned transient 2-7r?i/[/-periodic oscillations are vis- 
ible in the td dependence of the reflected field E rc n (t) at 
small t. 



VI. CONCLUSION 

In this paper wc generalized the familiar equilibrium 
expression for the optical conductivity in DMFT to 
the linear electromagnetic response of a nonequilibrium 
state. We find that the two-time optical conductivity 
cr(t, t'), which is probed in time-resolved optical spectro- 
scopics, can be expressed in terms of electronic real-time 
Green functions [see Eqs. (fT4|) and (fT8|) ]. which can be 
obtained from the DMFT solution. The expression for 
a(t, t') is completely general. Only anisotropic effects are 
disregarded that would lead to a dependence of the signal 
on the relative polarization direction of pump and probe 
pulses, i.e., averaging over the pump-probe direction is 
assumed. 

As a first application we have applied the theory to 
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a hypothetical pump-probe experiment on the Falicov- 
Kimball model. The pumping out of equilibrium was 
modelled by a sudden change in the interaction parame- 
ter, after which an electrical field pulse probes the relax- 
ation between metallic and insulating phases. We observe 
very fast relaxation with a relaxation time comparable to 
the inverse bandwidth, such that no clear separation of 
the time scales occurs between the intrinsic relaxation 
and electromagnetic response. Moreover, the two-time 
optical conductivity reveals transient oscillations in the 
response on a shorter time scale on the order of the in- 
verse interaction. These collapse-and-revival oscillations 
are expected to be very robust, e.g., for different densi- 
ties. Using time-resolved spectroscopy it may thus be 
possible to observe this phenomenon, which is known 
from experiments with ultracold atoms in optical lattices, 
in the relaxation of correlated electrons in solids as well. 

In the future, it should become feasible to solve the 
DMFT equations also for the Hubbard model in nonequi- 
librium. This will provide important insight into the 
dynamics of the pumped Mott insulator at short time- 
scales. 
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APPENDIX A: MOMENTUM SUMMATIONS 

For the homogeneous and isotropic relaxation without 
external fields discussed in Section [V] the evaluation of 
momentum sums is performed along the same lines as 
in equilibrium^ 6 - Because the DMFT self-energy Y, ka = 
S CT is local, the momentum k enters the DMFT equa- 
tions (Pl- tfTTj) only via the single-particle energy e ka a 
[Eq. (8b])], i.e., G kr7 {t,t') = G~ tk „ a (t,t') in zero fields 
The k sums in Eq. (JTT|, p8ap . and (|18b|) can then be 
reduced to integrals over a single energy variabl o 27 ' 50 by 
introducing the local density of states 



In Eq. (|A1[) . \ka) is the single particle state of the hop- 
ping matrix V?; for a Bravais lattice one has |(i|fccr)| 2 = 
1/N. For any function g(e) we thus obtain the relations 

oo 

= jde Pa (e)g(e) (A3) 

k -oo 

in Eq. (fTT|). and 

oo 

^E^-9(^)= J deD° p (e)g(e) (A4) 

k -oo 

OO 

jr^J2{dk a d kfj e ka )g(l kr7 ) = J de [d e D^(e)] g(e) 



(A5) 



in Eqs. (|18ap and (|18b[) . Here the last relation is proven 
using partial integration and the identity v ka d e 5(e — l ka ) 
= -d k 5(e - e k(J ). 

In this work we use a semielliptic density of states, 
Pl(e) = (2/-kW 2 )VW 2 — e 2 for the mobile particles in 
the Falicov-Kimball model, which leads to a simple 
self-consistency condition for the auxiliary single-site 
problem^ In the limit of infinite coordination number, 
this density of states is obtained for nearest-neighbor 
hopping on the Bethe lattice, but also for a particular 
choice of longer range hopping amplitudes on the hyper- 
cubic latticed In the latter case one obtains^ 



DUe) = <W 



Wa 2 



4/lVl - {e/wy 



cxp 



-2ci-r J 



ey/l - {e/W) 2 + WsisT 1 (e/W) 
nW/2 



(A6) 



for the dispersion function (|A2p . where a is the lattice 
constant. We adopt this form for the mobile particles 
in the Falicov-Kimball model; D a p = for the immobile 
species. 



(Al) 



and the dispersion function 
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